Anisotropic simple-cubic Ising lattice: 
extended phenomenological 
renormalization-group treatment 



M. A. Yurishchev * 

Institute of Problems of Chemical Physics of the Russian Academy 
of Sciences, 142432 Chernogolovka, Moscow Region, Russia 



Abstract 

Using transfer-matrix extended phenomenological renormalization-group 
methods [6] the improved estimates for the critical temperature of spin- 1/2 
Ising model on a simple-cubic lattice with partly anisotropic coupling strengths 
J = (J', J', J) are obtained. Universality of both fundamental critical expo- 
nents yt and yh is confirmed. We show also that the critical finite-size scaling 
amplitude ratios A^(i)An/A^, A^n /A^, and A^(a) /A^(A) are independent of the 
lattice anisotropy parameter A = J'/ J. 

PACS numbers: 05.50. +q, 05.70.Jk, 64.60.Fr, 75.10.Hk 



*Electronic address: yur@icp.ac.ru 



1 



1 Introduction 



The phenomenological renormalization-group (RG) method in which the transfer- 
matrix technique and finite-size scaling (FSS) ideas are combined is a powerful tool 
for investigation of critical properties in different two-dimensional systems [1, 2]. 
Unfortunately, its application in three and more dimensions is sharply retarded due 
to huge sizes of transfer matrices which arise in approximations of d-dimensional 
lattices by L^~^ x oo subsystems. 

Indeed, even in the simplest case of systems with only two states of a site (a spin- 

1 Ising model) the order of transfer matrix in three dimensions {d = 3) increases 
according to the law 2^^ (instead of the essentially more sparing law 2^ in two 
dimensions). Hence, for the cluster 3 x 3 x oo it is needed to solve the eigenproblem 
of the transfer matrix 512 x 512, for the 4 x 4 x oo subsystem — 65 536 x 65 536, 
and for the 5 x 5 x oo cluster it is required to find the eigenvalues and eigenvectors 
of dense matrices with huge sizes of 33 554432 by 33 554432. 

One can now solve the full eigenproblem for the transfer matrices of Ising par- 
allelepipeds L X L x oo with the side length L < 4. Our aim in this paper is to 
use such solutions with the largest effect and extract as much as possible accurate 
information about physical properties of the bulk system. 

The ordinary phenomenological RG is based on the FSS equations for correlation 
lengths [1, 2]. However, it is known [3, 4, 5] that the phenomenological RG can be 
built up using other quantities with a power divergence at the phase-transition point. 
It is remarkable that the such modified renormalizations can provide more precise 
results by the same sizes of subsystems [6]. 

In this article we give improved estimates for the critical temperature of the 
anisotropic three-dimensional (3D) Ising model by use of the extended phenomeno- 
logical RG schemes found in [6]. Achieved estimates are taken to calculate the values 
of different invariants of the 3D Ising universality class. Some of them are also given 
in the paper. 

2 Basic equations 

Start from the ordinary FSS equations [1, 2] for the inverse correlation length Kiit, h) 
and singular part of the dimensionless free-energy density /£(t, h), but we write them 
out for the derivatives with respect to the reduced temperature t = {T — Tc)/Tc and 
external field h: 

h) = u^y^+'^y^-^f/fit', h') (1) 

and 

/2^'"'")(i, h) = 6"^f*+"2//«-<i/2(^'")(i', h'). (2) 

Here KjJ^'^^\t, h) = d'"'^"-KL/ dt"'dh"' and the same for Z^''™'"''; Vt and Uh are, respec- 
tively, the thermal and magnetic critical exponents of the system; b — L/L' is the 
rescaling factor. Deriving Eqs. (1) and (2) we used a linearized form of RG equations 
f ~ b^H and h' ~ b^^h. 
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In the traditional phenomenological RG theory [1, 2], Eq. (1) with m = n = is 
considered as an RG mapping {t, h) {t', h') for a cluster pair (L, L'). By this, the 
critical temperature Tc is estimated from the equation 

Lkl{T,) = L'kl'{%). (3) 

Phenomenological renormalization (t, h) {f, h') can be also realized by using 
any of the relations (1) and (2) or their combination. It has been shown by the author 
[6] that some of such extended renormalizations lead to more rapid convergence 
in L than the standard phenomenological RG transformation. In particular, test 
examples on the fully isotropic systems [6] exhibited that the relations 



and 



Xl 



(4) ,,(1) 



(4) 



(5) 



L'^xl {L'Yxl 

locate Tc more accurately in comparison with the ordinary RG equation (3). In 
the relations (4) and (5), the derivative Kj^ = d'^Ki/dh'^, the zero-field susceptibility 
Xl = fi}'^''^\ and the nonlinear susceptibility Xl"* = /l*"^'^'' can be evaluated by 
standard formulas via the eigenvalues and eigenvectors of transfer matrices (see, 
e. g., [7, 8]). 

To get the thermal critical exponent yt we applied two approaches. Firstly, we 
used again the standard finite-size expression 

^ HLkL/{L'kL')] 

^* ln(L/L') ^ ' 

which follows from Eq. (1) hy m — 1, n — kl — duL/dt. Secondly, we took the 
formula 

^* " {KLKL'kLkL'Y/nniL/L'Y ^ ^ 

This expression is a direct sequence of the well-known Roomany-Wyld approximant 
to the Callan-Symanzik /3-function [2]. 

To calculate the magnetic critical exponent i/h we also used two ways: 

2HL/L') 

and 

^'^ 2^21n(L/L') 

[these finite-size relations follow from Eqs. (1) and (2)]. 

In addition, we calculated the universal ratios of critical ESS amplitudes. Such a 
kind of the ratios can be identified from the Privman-Eisher functional expressions 
[9] which for partly anisotropic systems read [8] 

KL{t, h) = L-^GolC{C\tLy\ C^hiy^) (10) 
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and 

fl(t, h) = L-'^GoJ^{CitLy\C2hLy^). (11) 

Scaling functions /C(xi,a;2) and T{xi^X2) arc the same within the hmits of a given 
universaUty class but they may depend on the boundary conditions and the sub- 
system shape (a cube, infinitely long parallelepipeds, etc.). All nonuniversality in- 
cluding the lattice anisotropy is absorbed in the geometry prefactor Gq and metric 
coefficients Ci and C2. The critical amplitude ratios from which the parameters Gq, 
Ci, and C2 drop out should be universal. In particular, the amplitude combinations 

(Binder-like ratio for the spatially anisotropic systems), 
and 

are expected to be not depend on the lattice anisotropy parameter A — J'/ J. 



3 Results and discussion 

In the present paper we carried out calculations for the subsystems L x L x 00 with 
L — 3 and 4. To avoid undesirable surface effects the periodic boundary conditions 
have been imposed in both transverse directions of parallelepipeds LxLx 00. Thus, 
the transfer matrices for which the eigenproblems was solved were dense matrices 
of sizes up to 65 536 x 65 536. To solve the eigenproblem we took into account the 
internal and lattice symmetries of subsystems and used the block-diagonalization 
method (see, e. g., [7]). Calculations were performed on an 800 MHz Pentium III 
PC running the FreeBSD operating system. 

Note that the matrix quasidiagonalization procedure is rather tedious in real- 
ization. Another possible way for treating the full eigenproblem of such large-scale 
matrices is to use supercomputers (say, Ref. [10] where the direct diagonalizations 
have been carried out for the Hamiltonian matrices of sizes 2^^ by 2^^). 



3.1 Critical temperature 

The critical temperature estimates coming from solutions of the transcendental equa- 
tions (4) and (5) are collected in Table 1. 

In the purely isotropic case ( J' = J) there are high precision numerical estimates 
for the critical point of the 3D Ising model. The most precise values for it have 
been obtained by Monte Carlo simulations [11, 12]; = 0.221654 59(10), i. e. 
kBTjJ = l/K^ = 4.511 5240(21). 
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Inspecting Tabic 1 one can see that the estimates for J' = J which follow from 
Eq. (4) and (5) are the lower and upper bounds respectively. By this, their mean 
value has the accuracy of 0.01%. Note also that our mean estimate is better than the 
value ksTc/ J — 4.533 71 obtained in Ref. [13] (see also [14]) for the fully isotropic 
lattice by using the ordinary phenomenological renormalization of the bars with 
L = 4 and 5. 

Discuss now the anisotropic case. Here there is well-known exact asymptotic 
formula for the critical temperature [15] 

{kBTc/J)asym = 2/[ln(J/2J') - lnln(J/2J') + 0(1)] (15) 

as J'/ J 0. It is a direct consequence of the molecular-field approximation in 
which the linear Ising chain is taken cluster. 

Unfortunately, simple formula (15) yields considerable errors in the region 10~^ < 
J'/ J < 1. Its modifications in spirit of Ref. [16], ksT^/J ~ 2/[ln(J/J')-lnln(J/J')], 
lead to a loss of monotonous convergence when J'/ J varies from unity to zero. 

We choose infinitely long chistcrs L x L x oo stretched in a lattice direction 
with the dominant interaction J. Such a cluster geometry refiects the physical 
situation in the system. Therefore one may expect more precise results for the 
critical temperature as the anisotropy of the quasi-one-dimensional lattice increases. 
We may also expect the monotonous convergence for the estimates from Eq. (4) 
and (5) because there are must be physical reasons (finite length of clusters in 
the longitudinal direction, etc.) for the non-monotonous or oscillatory character of 
behavior; they are absent in our approximations. That is, if Eq. (4) yields the lower 
bound in the most unfavorable case J' = J then it should preserve such behavior 
for all J' < J. Similar arguments are vahd for the estimates following from Eq. (5); 
these estimates are upper. 

Note that the mean values from Table 1 are better not only than the estimates of 
/cb^c/ J calculated with the (3, 4) cluster pair by the standard phenomenological RG 
method, but than their improvements found by means of three-point extrapolations 
from the sizes L — 2, 3, and 4 to the bulk hmit [17]. 

In the range 10~^ < J'/ J < 1, there are also the data for the critical temperature 
of a simple-cubic Ising lattice which were extracted from the Pade-approximant 
analysis of the high-temperature series [18]. For J' = J according to these data, 
ksTc/ J = 4.5106 that is lower by 0.014% in comparison with the results of Ref. [12]. 
For J' /J = 0.1 the authors of Ref. [18] found the value keTjJ = 1.343. This 
quantity overestimates somewhat the mean value from Table 1. At last, for J'/ J — 
0.01 the series method [18] yields ksTc/ J = 0.65 that goes out of our lower bound. 
This is not surprising because the calculations based on the high-temperature series 
rapidly deteriorate owing to the very limited number (< 11) of terms available in 
such series for the anisotropic lattices. 

So, we may treat the values found from Eqs. (4) and (5) as lower and upper 
bounds on the real critical temperature. Their mean value for each J'/ J yields the 
best estimate which we achieve in this paper for the reduced critical temperature 
ksTc/ J (the last column in Table 1). By this, its absolute error is not larger in any 
case than the half difference of the corresponding upper and lower bounds. Using 
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data from Table 1 we establish that the relative errors for kBTd J monotonically 
decrease from 0.72% to 0.14% as J' / J goes from 1 to 10~^. 

3.2 Invariants of the 3D Ising universality class 

Taking the improved estimates for the critical temperature of anisotropic simple- 
cubic lattice we calculate now some invariants of the three-dimensional Ising model 
universality class. 

3.2.1 Critical exponents 

According to the RG theory, critical exponents are determined entirely by a fixed 
point and do not depend on the lattice anisotropy. For a three-dimensional Ising 
model the universality of critical exponents has been confirmed for A G [0.2, 5] by 
the high-temperature series calculations [19]. 

At present, the most precise estimates of critical exponents are provided by the 
high-temperature expansions for models with improved potentials characterized by 
suppressed leading scaling corrections. For the fully isotropic 3D Ising lattice these 
methods yield v = 0.630 12(16) and 7 = 1.2373(2) [20]. Hence, = l/i/ = 1.5870(4) 
and yh = id + 7/z/)/2 = 2.481 80(18). 

In Table 2 we report our estimates for the critical exponents yt and y^- It follows 
from those data that when the lattice anisotropy parameter A varies in three orders 
(from unity to 10~^), the estimates of critical exponents are changed only on a few 
per cent or less. In particular, calculations by Eqs. (6) and (7) with the cluster pair 
(3,4) yield respectively yt = 1.47(6) and yt = 1.60(7). Their variations are over 
the range 4 — 4.4%. Similar calculations of the magnetic critical exponent carried 
out by use of Eqs. (8) and (9) also with the pair (3, 4) leads to yh — 2.586(5) and 
yh — 2.579(5), correspondingly. Relative dispersions of these estimates are about 
0.2%. 

Thus, our calculations confirm the universality of both critical exponents in a 
remarkably wider range of A than it was done in earlier investigations. Systematic 
errors of the achieved estimates arise due to small sizes, L, of subsystems used. 

3.2.2 Critical FSS amplitude ratios 

Critical amplitudes are determined by scaling functions. As a result, their "universal 
ratios" like A^(4)/^^(4) = /C^'^'^^(0, 0)/^^'^'^^(0, 0) depend, generally speaking, upon 
the lattice anisotropy because it can change the shape of subsystems. But in the 
case of parallelepipeds L'^^^ x 00 with unchanged (between themselves) transverse 
coupling constants the shape of a sample (all its aspect ratios) will be independent 
of the interaction in longitudinal direction. Such a kind of the universality is studied 
here. 

Table 3 contains results of our calculations for the critical FSS amplitude ratios 
U = A^(4,)Ak/A'^, Yi = A^"/A^, and Y2 = A^(4)/A^(4,). Calculations have been 
carried out for A G [10~^, 1] by use of a cyclic cluster 4 x 4 x 00. 
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In accord with data of Table 3, the average ratio U = 4.900(3). Hence, when 
the anisotropy parameter A varies in three orders, this quantity changes only on 
0.06%. With such accuracy we may consider the given ratio as a constant. In the 
case of fully isotropic lattice, A,^ — 1.26(5) and A^(4)/A'^ — 3.9(2) [8] and therefore 
A.^(4)Af,JA^ — 4.9(5). Our values of U from Table 3 are in good agreement with this 
estimate. 

It is follow from Table 3 that Yi = A^n/A^ = 1.759(2). Hence, the constancy of 
this universal amplitude ratio is estimated at least as a few times 10~^. Our average 
value for Yi agrees well with the estimate for the isotropic lattice, A^" /A^ — 1.749(6) 
[8]. 

According to data of Table 3 the amplitude ratio Y2 = A^{4)/A^(4) = 2.0133(6). 
Thus, this quantity is most stable out of all invariants of the 3D Ising universality 
class which were investigated in the paper. Note that we are not aware of any 
quantitative estimates for the Ai^(4.)/ A^(4). 

4 Conclusions 

In this paper the large-scale transfer-matrix computations of the phase-transition 
temperatures, as well as the critical exponents and critical FSS amplitudes have been 
carried out. Application of the extended phenomenological RG schemes allowed us 
to find the tight limits on the critical temperature in the anisotropic simple-cubic 
Ising lattice and improve the available estimates for it. 

We also calculated the thermal and magnetic critical exponents. Our results 
confirm the universality of yt within 4 — 4.4% and of yh within 0.2% over a wide 
range of A, 10"^ < A < 1. 

Finally, the presented results give an obvious evidence that the critical FSS 
amplitude ratios U — A^(4)Ai^/A^, Yi — A^n /A^, and Y2 — A^(a) j A^(a.) do not depend 
on the lattice anisotropy parameter A = J'/ J with accuracies at least 0.1%. We 
give likely for the first time in the literature an estimate for the universal quantity 
Y2. 
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Table 1: Lower and upper bounds on the critical temperature and their mean values 
(our estimates of ksTd J) in the 3D sc spin-| Ising lattice vs A = J' / J. Calculations 
with a cluster pair (3, 4). 



A Eq.(4) Eq.(5) mean 

LO 4.47965814 4.54424309 4.51195062 

0.5 2.91008665 2.94295713 2.92652189 

0.1 1.33649605 1.34570054 1.34109829 

0.05 1.03544938 1.04144927 1.03844933 

0.01 0.65054054 0.65323146 0.65188600 

0.005 0.55440490 0.55643112 0.55541801 

0.001 0.40743000 0.40859011 0.40801006 
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Table 2: Estimates of the thermal and magnetic critical exponents by different values 
of J' / J. Calculations with a cluster pair (3,4). 



Vt Vh 



J' /J 




Eq.(6) 


Eq.(7) 


Eq.(8) 


Eq.(9) 


1.0 


4.51195062 


1.5760695 


1.7246286 


2.5971647 


2.5886128 


0.5 


2.92652189 


1.5256373 


1.6636718 


2.5902006 


2.5819462 


0.1 


1.34109829 


1.4700811 


1.5972576 


2.5843305 


2.5766511 


0.05 


1.03844933 


1.4533899 


1.5791439 


2.5836720 


2.5761101 


0.01 


0.65188600 


1.4236178 


1.5480583 


2.5832982 


2.5758028 


0.005 


0.55541801 


1.4141719 


1.5383503 


2.5832029 


2.5757888 


0.001 


0.40801006 


1.3984754 


1.5222765 


2.5834573 


2.5757953 






1.47(6) 


1.60(7) 


2.586(5) 


2.579(5) 
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Table 3: Estimates of the universal critical FSS amplitude ratios U = A^(4)A^/A'^, 
Yi = A^(2)/A^, and Y2 = A^(4)/A^(4} for the Ising system with the cylindrical geom- 
etry L X L X 00 and periodic boundary conditions. Data for L = 4. 



J'/J 


ksTc/J 


U 




Y2 


1.0 


4.51195062 


4.8956599 


1.7550004 


2.0146443 


0.5 


2.92652189 


4.8967625 


1.7572512 


2.0136519 


0.1 


1.34109829 


4.9011909 


1.7596003 


2.0129829 


0.05 


1.03844933 


4.9014406 


1.7597697 


2.0129285 


0.01 


0.65188600 


4.9015375 


1.7598563 


2.0128977 


0.005 


0.55541801 


4.9015529 


1.7598646 


2.0128953 


0.001 


0.40801006 


4.9015782 


1.7598732 


2.0128938 






4.900(3) 


1.759(2) 


2.0133(6) 
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